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Abstract 


A best fit ellipse and hyperbola have been calculated to represent 
several hundred magnetopause and bow shock positions observed by six IMP 
spacecraft* Average geocentric distances to the magnetopause and bow 
shock near the ecliptic plane are 11.0 R E and l4 .6 Rg in the sunward 
direction ? 15.1 R e and 22.8 R E in the dawn meridian and 15.8 and 27.6 R E 
in the dusk meridian. The bow shock hyperbola is oriented in a direction 
consistent with that expected considering aberration of a radial solar 
wind. Observed magnetopause crossings agree well with theoretical 
predictions in the noon meridian plane but fall outside the theoretical 
boundaries in the dawn -dusk meridian planes. IMP 4 plasma data are used 
to demonstrate that the solar wind momentum flux is the prime factor 
controlling the orbit-to-orbit changes in the boundary positions. Data 
suggest that the interplanetary fieH orientation also affects the 
distance to the magnetopause boundary with more earthward crossings 
corresponding to southward fields. Six unusual bow shock locations up 
to 22 R e beyond the average position are found to be due to an enhanced 
standoff distance associated with a low Alfven Mach number. The possibility 
is raised that the solar wind may have become sub-Alfvenic on July 31* 

1967. 
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Introduction 

The earth’s magnetopause and how shock have been detected by numerous 
types of experiments on virtually every scientific spacecraft traversing the 
appropriate regions of space. The spacecraft Pioneer 1 (Sonett et al., i960), 
Pioneer 5 (Coleman, 1964), Explorer 10 (Heppner et al,, 1963 5 Bonetti et al. , 
1963) and early Soviet probes (Gringauz et al., I96I; Shklovskii et al., 
i960) provided the early measurements in the boundary regions but the 
exploratory nature of the experiments plus limited quantities of data prevented 
a clear understanding of the underlying physics of the solar wind-earth 
interaction. The first definitive studies which included mapping of the 
positions of the boundaries were carried out with data from the spacecraft 
Explorer 12 (Chaill and Amazeen, 1963; Freeman, 1964; Cahill and Patel, 1967) 
(magnetopause only) and IMP 1 (Ness et al., 1964; Bridge et al., 1965; 

Wolfe et al., 1966) (magnetopause and shock). Subsequently the spacecraft 
Explorer 14 (Frank and Van Allen, 1964), IMP 2 (Fairfield and Ness, 196?; 
Binsack, 1968) IMP 3 (Ness, 1967), OGO-1 (Heppner et al., 1967; Holzer 
et al., 1966), Vela 2A and 2B (Gosling et al., 1967) and 0G0-3 (Russell et 

’ ‘ i 

al., 1968) have further refined these measurements within 35 Re and the 
spacecraft Explorers 33 and 35 (Behannon, .1968, I97O5 Mihalov et al., 1970; 
Howe, 1970) have exten^sd the observations to greater distances behind the 
earth. 

The general picture revealed by all these measurements is of time , 

dependent magnetopause and shock positions whose distances from the center 
of the earth are scattered about average values which are near 11 R E and 
l4 R e respectively in the solar direction and 15 R E and 25 R £ in the meridian 


planes. To a first approximation the boundaries are symmetrical about 
the earth-sun line but due to the earths motion about the sun a 3°-5° 
deviation from symmetry is expected. Walters (1964) predicted an additional 
asymmetry due to the interplanetary field wh-ch on the average should increase 
the asymmetry due to the orbital motion. Those spacecraft making both dawn 

and dusk hemisphere measurements (Heppner et al., 1967; Ness, 1967? Gosling 

} 

et al., 1967; Behannon, I968.J Mihalov et al., 1970) confirm an asymmetry but 
it is not necessarily larger, than that expected for aberration. A larger 
8° skewing suggested by Hundhausen et al. (1969) on the basis of observed 
flow directions -in the magnetosheath has not been confirmed by boundary 
position measurements. 

The magnetopause crossing nearest the earth is that detected by the 
ATS 1 spacecraft at 6.6 Rg (Opp, 1968, and companion papers) which was 

accompanied by an abnormally close in bow shock at 13.4 Rg near the dawn 

» 

meridian (Russell et al., 1968). An abnormally distant bow shock location 
(Heppner et al., 1967) occurred approximately 4 Rg beyond the average 
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position and was explained by magnetosphere inflation plus an unusually low 
Alfven Mach number which theory (Spreiter and Jones, 1964) predicts 
should enhance the standoff distance of the bow shock. Gosling et al. (1967) 
have presented several instances when observed plasma parameters changed 
in a manner that was consistent with boundary motion and world wide geomagnetic 
field compression events. Binsack and Vasyliunas (1968) had rather good 
success predicting the observed position of the boundaries using solar wind 
data obtained at the time of a boundary crossing. This suggests that boundary 
positions are controlled primarily by the solar wind momentum flux which 
compresses the magnetosphere to a greater or lesser extent at different 
times. 
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Further studies of solar 'wind control of boundary position have 
generally utilized the reported positive correlation of solar wind 
velocity and geomagnetic activity index Kp (Snyder et al., 1963). Since the 
solar wind velocity tends to be high when Kp is high, several workers 
(Patel and Dessler, 1966; Cahill and Patel, 1967$ Holzer et al,, 19665 
Heppner et al,, 1967? Gosling et al., 1967) have searched for a correlation 
between boundary position and Kp. The net result of these studies is that 
there is a weak correlation between boundary position and Kp with a 
tendency for more earthward positions to correspond to high Kp. Lack of a 
better correlation is probably due to the fact that there exists an inverse 
relationship between solar wind density and velocity (Neugebauer and Snyder, 19665 
Hundhausen et al., 1970 5 Burlaga and Ogilvie, 1970a) which tends to keep the 
flux constant even though the velocity changes. Snyder et al., (1963) and 

Neugebauer and Snyder (1966), in fact* report that momentum flux does not 
correlate as well with Kp as does velocity. Recently Meng (1970) has 
demonstrated that distant magnetopause crossings invariably correspond to 
quiet conditions (low AE index), whereas crossings nearer the earth may 
correspond to quiet or disturbed conditions. Aubry et al. (1970) report a 

case where the solar wind flux remains constant but the magnetopause moves 

* 

inward as the interplanetary field becomes southward. Their suggestion 
is that the interplanetary field direction exerts some control on the boundary 
position by eroding magnetic flux from the subsolar magnetosphere. 


Analysis 


The present study was undertaken to establish an accurate representation 
of the average position and shape of the magnetopause and bow shock and to 
investigate the magnitude of the variations from the average and. their 
causes. For this purpose, plots of magnetic field direction and magnitude 
versus time from the six spacecraft IMP*s 1-4, Explorer 33 and Explorer 35 
were scanned and the magnetopause and bow shock positions were tabulated for 
each pass. The distribution of 3^9 shock crossings and 4?4 magnetopause 
crossings by spacecraft is given in Table I. The positions determined in cases 
of multiple boundary crossings were average positions in the sense that the 
interval of magnetosheath (magnetosphere) data on the sunside of the selected 
bow shock (magnetopause) location was chosen to be e<£ual to the interval of 
interplanetary ^magnetosheath) data inside the selected bow shock 
(magnetopause) location. Boundary positions were listed in the solar 
ecliptic coordinate system. This system was found to be approximately as 
good as solar magnetospheric coordinates and better than solar magnetic 
coordinates in ordering Vela data (Gosling et al., I967). 

Shock crossings are characterized by field magnitude increases by a 
factor which is typically 2-4 and occurs over a time interval which is 
usually small compared to adjacent plotted points, (see Table I for the number 
of points measured and plotted each minute). The bow shock was observed on 
every pass when data was available in the appropriate region although occasionally 
the exact position was uncertain by a few tens of minutes (a few tenths of 
an earth radii) due to rapid multiple crossings, upstream waves or unusual 
shock thicknesses. Orbits where the apogee was not well beyond the average 
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shock position were omitted (all of IMP 2) so as not to bias the sample 
toward the more earthward crossings. Intervals when the unique Explorer 33 
trajectory ran approximately parallel to a boundary were also omitted. 

Magnetopause crossings are characterized by discontinuities in field 
direction and magnitude (e.g., Hyde, 1967) and a change in the field 
fluctuation level. Boundary identification occasionally becomes difficult 
when the magnetosphere and magnetosheath fields happen to align themselves 
and on approxina tely 5$ of the passes a boundary could not be chosen with 
reasonable confidence. This problem was somewhat more prevalent on the 
earlier spacecraft with lower sampling rates. Explorer 35 crossings of 
the magnetic tail boundary have been reported by Mihalov et al. (1970) and 
were not reinvestigated. 

In order to obtain an accurate analytical representation for the bow 
shock and magnetopause a computer program was written to obtain the best fit 
conic to the two-dimensional representation of the magnetopause or shock 
data in the solar* ecliptic plane. This program in effect translates, rotates 
and changes the shape of a conic in order to minimize the sum of the 
differences between the data points and the curve. Only the type of conic 
was specified for a given run (an ellipse for the magnetopause, a hyperbola 
for the bow shock) along with reasonable starting values which were found not 
to effect the final solution. For computational ease a convex body norm 
(Householder, 1958) was used whereby the distance between data points 
and fitted curve was calculated along the line between the data point 
and the point (-5,0) or (-30,0) for the ellipse and hyperbola respectively. 

By imposing no constraints on the center or orientation of the conic 
it was hoped that an east-west asymmetry could be detected. 
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Although there is no guarantee c? a unique solution to such a curve fitting 
problem, the data was well enough ordered such that the program 
invariably approached a solution that was reasonable and the best fit 
for a limited region of parameter space. Results are given in Section 3, 
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Results 

Average Boundary Positions 

The crosses in Figure 1 represent the average position of the bow 
shock on 188 passes with position j Z se |<7 R E . Figure 2 shows similar 
data for 255 passes through the magnetopause region. The actual three 
dimensional locations have been converted to two dimensions by rotating 
the points into the solar ecliptic plane. The rotation is about the 
X ge axis for X ge <0 and in a meridian plane for X xe >0, This means that 
the two dimensional representation is valid only to the extent' that the 
boundaries exhibit spherical and cylindrical symmetry in the subsolar and 
antisolar hemispheres respectively. Each poir«‘ in the ecliptic plane 
is then converted to the prime coordinate system by a rotation of 4° 
to eliminate the expected aberration due to the earth* s motion in the 
presence of an average 420 km/sec solar wind. 

The best fit conics (see Section 2) are illustrated by the solid 
curves in the figures and can be expressed by the equation 

2 2 

y + Axy + Bx + Cy + Dx .+ E - 0 (1) 

with the constants given in Columns 1 and 4 of Table, 2. The tail boundary 
crossings of Mihalov et al. (1970) have been rotated by 4 ° and plotted 
in Figure 2 but were not used in the curve fitting. Straight lines have been 
joined to the ellipse at X = -15 R E to represent the extended tail. The 
dashed shock curve in Figure 1 is the theoretical curve obtained by 
Spreiter and Jones (1963) using the dashed theoretical magnetopause 
of Beard (i960) . The dashed magnetopause in Figure 2 represents the 
more recent work of Olson (I969). The geocentric distances to the 
experimental curves at the sub solar point are 10. 9 and l4.5 for the 
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magnetopause and shock respectively. The magnetopause (shock) 
intersections with the dawn and dusk meridian planes are -15*7 (-24.9) 
and 15.3 (26.2) respectively. Agreement between experiment and theory 
is quite good with the primary discrepancy being the increased dimension 
of the experimental magnetosphere in the dawn dusk plane. The 
experimental deviation from symmetry after the 4° rotation is apparent 
from the difference between the solid curve and the symmetrical 
theoretical dashed curve. 

Best fit curves to the data analyzed in a different manner have 

the constants listed in Columns 2-3 and 5-6 of Table, 2. Columns 2 and 5 

describe curves fit to data which was rotated to the meridian plane 

as described above but for which the 4° aberration correction was 

not performed. Column 3 and 5 represent curves fit to data which has all 

been rotated into the ecliptic plane about the X axis and for which 
o 

no 4 rotation has been performed. The fourth row of the table lists 
the angle between the X axis and the axis of the conic. The angle is 
defined as positive in the direction east of the sun or in the sense 
of an increasing aberration effect. Columns 2 and 3 show the 

^ o « 

expected 4 difference from column 1 and the absolute values suggest 

o , " ' : . 

an angle approximately 1 larger than expected from aberration. The 
orientation of the magetopause' fails to show any rotation in the direction 

expected due to the earth’s motion about the sun. The net effect is that in the 

. o 

4 aberration corrected system the magnetosheath -in the dawn-dusk 
meridian plane is more than 1.5 % or w i fier in the dusk quadrant 
than in the dawn quadrant. In an uncorrected coordinate system where the 
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dawn dusk meridian plane cuts asymetrically through the boundaries 
this difference is about 4 R E or almost 50 1o wider in the dusk meridian. 

Since dowrurtream boundary measurements at the moon’s orbit (Mihalov et 
al., 1970) and measurements within the magnetotail (Behannon, 1970) 
clearly revealed the abberation effect, either the magnetopause 
measurements nearer the subsolar hemisphere are not sensitive enough 
to reveal aberration, or else asymmetries within the magnetosphere 
are contributing an effect. Heppner et al., (1967) have made the latter 
suggestion in proposing that n bumps ft may form on the magnetopause 
in response to variations of the magnetosphere plasma pressure with 
latitude and longitude. 

Another fact to be deduced from Figures 1 and 2 concerns the 
effective ratio of specific heats, 7. Spreiter et al. (1966) have 
noted that the ratio of the shock standoff distance (subsolar shock 
distance minus subsolar magnetopause distance) to the subsolar 
magnetopause distance is essentially independent of Mach number 
for large Mach numbers and dependent only on 7. Consequently, if the 
solar wind fulfills the large Mach number criterion the experimental value 
of the standoff distance can be used as an indicator of the appropriate 7. 

The frequency distribution of Mach numbers is presented in Figure 3 
which was prepared using the IMP 4 plasma data of Ogilive and Burlaga (1970) 
along with magnetic field data. All available hourly averages of plasma 
density, velocity, proton temperature and field strength have been used to 
calculate both the Alfven Mach number M A « V/V A where V is the solar wind 
velocity and V A is the Alfven velocity (H^/4rrp) 2 where H is the measured 

' ’ ' ' ■ ■ . . , . /■ . ■ ‘ ' r i 

magnetic field strength and p is the mass density computed from the observed 
density of protons. The gasdynamic Mach number is M = v/a where 


where p 0 and p p represent the proton and electron pressures, p is the 
mass density, k is the Boltzman constant, nip is the proton mass, and T^ 
is the proton temperature. The electron temperature T is not 
measured on IMP 4 but it has been shown (Montgomery et al. , 1968; Burlaga 
and Ogilvie, 1970b) to be approximately independent of solar wind 
conditions at a value of 1.5 x 10^°K which is the value used here. It is 
not completely clear which is the appropriate Mach number but Spreiter 
et al. (I 966 ) suggest that M is appropriate when the magnetic field 
alignment is arbitrary and the field magnitude is small enough such that 
>>1. Figure 3 indicates that both M and M^ are usually large 
enough such that the average experimental value for the standoff 
distance should be independent of Mach number to the extent that it 
can be used as an indicator of the appropriate 7 . 

Table 2 indicates that the ratio of standoff distance to magnetopause 
distance has a value 0.33* According to Figure 16 of Spreiter et al. 

(1966) this falls between the values expected for 7=5/3 and 7 *2 
though it is slightly nearer to the 7 = 2 value. It should be, remembered, 
however, that if occasional low Mach number shocks could be eliminated 
from the experimental data the standoff distance would be reduced and a 
slightly smaller 7 would be selected. Also if Spreiter et al. had used 
a wider magnetopause, as now appears appropriate, the standoff distances 
would be slightly increased, their curve would be raised and a smaller 7 
would be indicated. In spite of these restriction 7 f s less than 
5/3 appear to be inappropriate for use in the solar wind. 
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Magnetopause crossings within 15 ° of the noon -midnight meridian plane 
and the dawn -dusk meridian plane are shown in Figures 4 and 5 respectively. 
Solid lines represent theoretical boundaries (Olson, 1969) for three 
different values of the geomagnetic latitude of the subsolar point, 
Olson 1 s boundary was computed for a subsolar distance of 10,7 B E which is 
near the experimental value obtained in this paper. The agreement with 
theory is quite adequate in the noon-midnight meridian plane considering 
the wide range of solar wind conditions and dipole orientations occurring 
for the measured points. The positions within 15 ° of the dawn and dusk 
meridian plane have been superposed in Figure 4 . In this plane it is clear 
that the great majority of the experimental points are located outside the 
theoretical boundary. This discrepancy is probably due to the 

presence of plasma in the magnetosphere (Heppner et al., 1967; Vasyliunas, 
1968) and the presence of the bow shock, both of which are completely 
neglected in the theory. 
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Normal Boundary Variations 

In order to test the hypothesis that a variable solar wind is 
responsible for the differences in boundary position, the IMP 4 interplanetary 
plasma measurements were analyzed in conjunction with the observed boundary 
crossing. Theory predicts that the subsolar magnetopause distance is given 
by the formula 



(e.g., Spreiter et al., 1968) where H 0 = .312 Gauss is the geomagnetic field 
strength at the earth's surface on the equator and p, the mass density, 
is here taken as plasma density n times proton mass mp. The factor f 2 /K 
is deteimtined by the physics of the interaction (Spreiter et al., 1 966; 
Schield, 1969). The factor K in the denominator is a measure of how 
efficiently solar wind particles transfer their momentum to the magnetosphe’ 
(K = 1 for inelastic collisions, K = 2 for elastic collisions, K » .8 in 
gasdynaiuic theory). The factor f relates the geomagnetic field Just inside 
the magnetopause to the undistorted dipole field at that point. The 
quanity f probably assumes a value between 1 (image dipole model) and 1.54 
(in a model with a ring current (Shield, 1969)). The factor f /K is taken 
to be unity in the present calculation. 

Using the IMP 4 interplanetary plasma data and equation 2 a subsolar 
magnetopause distance was obtained for each hour the spacecraft was in the 
interplanetary medium. In addition a prediction of the shock position was 
obtained using the work of Spreiter and Jones (1963). Figure 1 in their 
paper predicts that the shock standoff distance should be approximately 1.37D 
for solar wind mach numbers greater than 10 and successively, higher for 
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decreasing Mach numbers. The gas dynamic Mach number M rather than 
Alfven Mach number has been used in the present work in accord with a 
later suggestion (Spreiter et al., 1966 ). 

The magnetopause position predictions are illustrated by the lower 
trace on each grid in Figure 6 , Shown above the magnetopause prediction 
is the shock prediction given by I. 37 D. Nearly coincident but just above 
this trace is the more refined Mach number -dependent prediction. The close 
proximity of these upper two traces reflect the fact that M is generally 
larger than 5 (see Figure 3). Also shown in Figure 6 are the experimentally 
observed positions of the bow shock (circles) and magnetopause (crosses). 

In order to normalize the observed position data and eliminate the effect of the 
flaring out of the boundaries with longitude, differences between the observed 
boundary positions and the average fitted curves are plotted relative to the 
horizontal lines representing the average subsolar distances. The data shown 
is all from within approximately 50° of the earth-sun -line. Geomagnetic 
conditions were relatively quiet during this interval with four small SC 
storms occurring (marked by triangles in Figure 6 ) of which only one 
(September 19-21 ) had a significant Dst which was not larger than 80y 
(Sugiura and Cain, 1970). It is apparent in Figure 6 that the bow shock 
observations agree quite well with the predictions but the magnetopause 
observations are frequently substantially different from the predictions. 

Since the shock predictions are highly dependent on the magnetopause 
prediction it can be Concluded that the solar wind data accurately predicts 
the position of the average magnetopause. The greater variation in observed 
and predicted magnetopause positions relative to the observed and predicted 
shock positions may be due to (l) the fact that the solar wind may change 
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in the several hours between the last interplanetary measurement and the obser- 
vation of the magnetopause, (2) irregularities on the surface of the 
magnetopause producing variations in the 3.onally observed boundary position 
which do not substantially affect the shock location, and/or (3) the 
fact that the observed magnetopause position is more likely to be inaccurately 
identified since it is more difficult to determine than the bow shock. Since 
possibilities one and three could be eliminated with simultaneous interplanetary 
and boundary data and multi -experiment detection of the magnetopause, the 
above technique offers a possible means of detecting the presence of local 
irregularities on the magnetopause in future work. From Figure 6 it can be 
concluded that the position of the bow shock can be predicted to better than 
1 % 8o$ of the time and to better than .5 % 50$ of the time. 

Variations in the observed positions of the magnetopause by all the 
spacecraft are compared to the variations predicted from the interplanetary 
plasma data in Figure 7* The distribution of variations of 137 observed 
magnetopause positions from the fitted curve is represented by the dashed 
line and the prediction of the IMP 4 plasma da.ta are represented by the 
solid line. Zero variation from the fitted curve has been labeled 10.9 R„. 

ht 

The plasma predictions have been centered on this value even though the 
average distance predicted by the plasma data with f 2 /K = 1 was 10.3 Rg. 

The difference could be due either to the improper f^/K or a systematic 
error in the density which could be as large as + ,30$ (Ogilvie et al., 

1967). To bring the predicted values into agreement with the observations 
a 30$ decrease in density is needed, however there are other indications 
(Burlaga and Ogilvie $ 1970b) that the IMP 4 densities should not be 
reduced. It appears therefore that the likely source of the discrepancy is 
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in the value f 2 /K. Accepting the measured average values of n and v 
requires an f 2 /K of 1.4 to bring the predictions into line with the 
observations. The somewhat broader width of the observed distribution in 
Figure 7 again suggests that effects other than a uniform pressure balance 
may control the position of the magnetopause. It is also possible that 
a solar cycle variation in the solar wind momentum flux is broadening the 
distribution of experimental points. It is very difficult to test for 
such a variation because the measurements from each spacecraft 
( each year) are made at characteristic latitudes and with certain dipole 
inclinations at certain longitudes. All this results in a complex combination 
of several small effects which may influence boundary position. 

To pursue the relation between boundary position and magnetic activity 
the analysis of Meng (1970) was repeated using more than three times as much 
data as was previously available. A plot of the geomagnetic AE index 
(Davis and Sugiura, 1966) vs. the magnetopause position relative to its 
average location confirmed the result of Meng that an abnormally earthward 
magnetopause may correspond to either quiet or disturbed geomagnetic 
conditions but a distant magnetopause invariably corresponds to quiet periods. 
Because of the great similarity to Meng's Figure 1 this data is not presented 
here . 

In order to pursue the suggestion of Aubry et al. (1970) that 
interplanetary orientation affects magnetopause location, all 178 of the 
available IMP 4 magnetopause crossings were utilized in a study of boundary 
position as a function of field orientation outside the boundary. For 6 % 
of the crossings it was judged that the magnetosheath field during the 
approximately 2 hour period the spacecraft was nearest the magnetopause 
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was fairly steady and could be characterized as being predominantly northward 
or southward. Two groups of northward ( 62 ) and southward ( 54 ) crossings were 
thus determined and ellipses were fit to these two groups. The subsolar 
distances were 10.5 +1.3 and. 11.6 +2.0 for the south and north groups 
respectively where the plus or minus figure represents the standard 
deviation of the points from the fitted curve # Since the variability is 
due largely to solar wind pressure changes (f this 1.1 Rg difference in 
the average position would seem quite significant. When the analysis was 
repeated with IMP 1 > 2 and 3 data a similar result was obtained although 
the difference between the subsolar points was only 0.3 R^, This relation 
between field direction and boundary position supports the suggestion of 
Aubry et al. , that the magnetopause moves inward when the interplanetary 
field is southward. 

It is well-known that geomagnetic disturbance is associated with a 
southward directed interplanetary magnetic field (e.g., Hirshberg and 
Colburn, 1969) as well as a high solar wind velocity. Apparently the 
weak tendency for earthward boundary crossings to be associated with 
high Kp is also due to a southward field and not only to enhanced momentum 
caused by higher velocities. 


4 
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Distant Bow Shock Observations 

In the course of determining bow shock locations, six passes were discovered 
where the bow shock was found to be at least 10 R^ beyond the average location. 
These locations are designated by line segments in the upstream region 
in Figure 1 and they are labeled with the number of crossings witch occurred 
during the interval. These events are listed in Table 3 which includes the 
observing spacecraft, date and time of the first crossing and duration of 
the interval of crossings, number of crossings, average solar ecliptic position, 
average Kp and Dst values for the interval and average of the interplanetary 
field magnitude observed outside the various shock crossings. Clearly Kp 
and Dst bear no important relation to the abnormal locations. In the magneto- 
sheath adjacent to these crossings the field tended to be very quiet 
(Fairfield and Ness, 1970). 

In the first three of the IMP 4 distant shock events listed in Table 3> 

IMP 4 plasma data is available in the upstream region adjacent to the shocks 
and it may be used in trying to explain their distant locations. For two of 
these passes MIT Explorer 35 interplanetary plasma data was generously 
supplied by J. Binsaek and H. Howe for the entire interval of the pass-. 

During the July 5 event the IMP 4 plasma experiment measured plasma 
densities between 1.0 and 1.5 particles/cc and solar wind velocities between 
420 and 440 km/sec in the interplanetary region adjacent to the distance 
shock crossings. These parameters predict a subsolar magnetopause distance 

• . ; ' • ■ • " fj. 
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near 12.5 R E (Equation 2) and subsolar shock distance (1.37D) of approximately 
17 R e if the Mach number i»* high. The corresponding distance to the shock 
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at the subsatellite point is 25 which is 6-7 R E inside of the observed 
position of the shock. The Airven Mach number at the time of these distant 
shock crossings is between 1.2 and 1.6 which according to the work of 
Spreiter and Jones would account for the observed increase. A small ring 
current (Dst « -10y) would also tend to increase the boundary distances 
by a small factor. 

On July 30 at 23:20 the IMP 4 spacecraft observed the bow shock at a 
position 17 Eg beyond the average position of the shock. Both IMP 4 and 
Explorer 35 plasma experiments measure a solar wind with velocity of 325 km/sec 
and density .5 particles/cc at the time of the shock crossing. These values 
predict a subsolar magnetopause distance of 16.1 K E and a subsolar shock 
(1.37P) distance of 22 R E . This subsolar distance corresponds to a 
subsatellite distance of 23 R E which is at least 10 R E inside the observed 
location. With this plasma density of .5 and the observed interplanetary 
field magnitude of 7.5y the Alfven mach number is 1.4 which again is 
unusually low and can easily explain the large distance of the observation. 

During the 24 hours following this shock crossing the inbound IMP 4 
spacecraft observes a magnetic field which remains enhanced relative to 
the interplanetary field measured further upstream (X qe - 23, Y eo = -51, Z =3) 
by Explorer 35 by a factor 1.4 to 1.7. In the interval between 01:00 and 
23:00 on July 31? the plasma flux at Explorer 35 reached an even lower 
level that in fact was below the detectability level of the experiment. 

It is interesting to note that if this decrease to an unobservable level 
was caused by a density decrease of a factor of 2 from the density value 
of .5? or a 307. decrease in velocity, the solar wind becomes sub-Alfvenic. 
Throughout this interval of possible sub-Alfvenic solar wind the IMP 4 
magnetic field strength 
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remained very quiet and greater than that measured further upstream 
"by Explorer 35 by a factor of 1.4 to 1.7. The enhanced magnitude at 
IMP 4 could indicate the continued presence of an upstream shock but it 
could also be the normal increase associated with a new shockless mode 
of interaction between the geomagnetic field and a sub-Alfvenic 
solar wind. 

The abnormal shock locations on November 27 are somewhat more difficult 
to explain. Measured plasma parameter by Explorer 35 and IMP 4 at the time 
of the initial and innermost shock crossing indicate that shock location 
should be very near the average position if the Mach number is high. The 
Alfven Mach number is in fact 3 which, although unusually low, should increase 
its distance by only about 1 R E . The observed position remains about 5 Rjj. 
beyond this predicted location. The solar wind is measured on Explorer 35 
to be coming from 10° north of the ecliptic plane at this time which 
should rotate the shock outward at the position of the spacecraft which is 
6.5 Re below the ecliptic plane. A small ring current would also tend to 
increase the observed location. On subsequent crossings during this pass 
at greater distances the density observed at IMP 4 decreases as low as 1.5 
parti cles/cm^ so that the predicted shock moves outward about 3 Rg» She 
magnetic field decreases to about 6y, however, so that the Alfven Mach 
number remains at a value near 3 . The discrepancy between prediction and 
observation remains as large as 8 Rg unless the secondary effects are considered. 

Although no plasma data are available for the remaining three events 
in Table 3* it should be noted that each of these events has an unusually 
strong interplanetary field strength associated with it. Since Alfven 
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Mach number ic inversely proportional to the field strength, these events 
are consistent with the suggestion that an abnormally low Alfven Mach number is 
responsible for an increased standoff distance and the unusually distant 
location of the bow shock. 


The frequency distribution of Alfven Mach numbers shown in Figure 3 


illustrates how unusual low Mach numbers are. It- nan be seen that Alfven 
Mach numbers less than 3 constitute less than 2.2$ of the total and those 
less than 2 (the July 5 and July 30 distant shock events) are hardly ever 
seen by earth orbiting spacecraft. Field strengths greater than lOy (the 
3 events without plasma data) are observed less than 10$ of the time and 
those greater than 15y (November 4, May 3 events) less than 2$ of the time 
(Ness, 1969 ). 

To obtain evidence that shocks located at greater distances from the 
earth are unusually weak shocks associated with lower Mach numbers the increase 
in field strength across the shock was investigated. From the MHD shock 
equations it can be shown (Whang, private comnunication) that the ratio 
of the jump in tangential field components across a perpendicular shock 


is given by the expression 



(2 1^+ 5P+5) + 


(2 ^ + 5P + 5) 2 + 32 M 2 ] 

16 4 


1/2 


(3) 


where and Hg are the upstream and downstream tangential field components 


respectively and £ is the ratio of P^ 4* to 


Id energy H / 


8tt. This 


ratio is relatively insensitive to £ and is a decreasing function of 

with the limit 1/ -*1/4 as M -* «. Although the actual shocks are not f 

Hg A 

necessarily the perpendicular shocks to which the equations refer, this 
decreasing tendency should carry over to the more general case. To compare 
the data to the predictions of this equation a simple hyperboloid of 
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revolution was used to locate the approximate plane of the shock and the 
measured field components parallel to this plane were determined. For 
the distant shock crossings the tangential field component increases were 
found to range from 1.4 to 2.5. These numbers correspond to Mach numbers 
below 3,3 according to equation 3. For 77 shocks at more normal positions 
throughout the subsolar hemisphere the corresponding field component 
Increase was by a factor of 3.2 which corresponds to a Mach number of 5. 
This result supports the conclusion that the distant shocks are weak and 
associated with low M&ch numbers. 


Summary and Conclusions 


Magnetopause and bow shock locations measured by six IMP spacecraft 
have been analyzed to determine* 1 the average boundary positions and the 
causes for their variations with time. Best fit curves obtained to 
represent the average boundaries in the solar ecliptic plane are 
characterized by geocentric distances to the magnetopause and bow 
shock of 11.0 and 14.6 R E respectively near the subsolar point. The 
average bow shock orientation is symmetrical about the expected incident 
direction of the solar wind to better than 2°. The average magnetopause 
orientation deduced from tubsolar hemisphere measurements, on the other 
hand, fails to reveal a corresponding aberration effect though such an 
effect is clear from tail measurements in the downstream region. The 
width of the average magnetosheath in the solar ecliptic dawn dusk plane 
is almost 50 $ greater in the dusk hemisphere as compared to the dawn 
hemisphere because of asymmetry of the shock. The usual theoretical - 
methods of calculating the shape of the magnetopause appear to be 
adequate in the noon-meridian plane but predict a magnetopause which is 
too close to the earth in the dawn-dusk meridian plane. This discrepancy 
is probably due to the theoretical assumption that there is no 
bow shock and no plasma in the magnetosphere. It is suggested that the 
experimentally determined value of .33 for the ratio of the shock standoff 
distance to tl)e geocentric subsolar magnetopause distance along with 
relatively high values of the gasdynamic Mach number together imply that an 
effective value of y between 5/3 and 2 is appropriate for the solar 
wind- earth interaction. 
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Analysis of IMP 4 plasma data in conjunction with the boundary 
positions shows that solar wind flux appears to be the primary factor 
controlling the average position of the magnetopause and bow shock observed 
on any given spacecraft pass. Knowledge of solar wind density and velocity 
is found to be adequate to predict the average position of the bow shock 
to better than 1 Rg 8 of> of the time and 0.5 Hg 50 io of the time when the 
geomagnetic ring current is not unusually large. Evidence is presented 
suggesting that the factor f^/K used in predicting the magnetopause position 
assumes a value greater than unity. Prediction of the exact magnetopause 
position may be more difficult than prediction of the shock position if 
local variations of the boundary are important. Analysis confirms Meng's 
(1970) finding that distant magnetopause positions correspond to quiet 
condition^ whereas earthward positions are observed during either disturbed 
or quiet times. The direction of the interplanetary magnetic field is 
found to be a secondary factor influencing the boundary positions with a 
southward field producing a more earthward magnetopause. This supports 
the suggestion of Aubry et al. (1970) that in the presence of a southward 
interplanetary field magnetic flux is eroded from the subsolar magnetosphere 
and added to the tail. 

On rare occasions the bow shock is observed at exceptionally distant 
locations, as much as 22 Rg beyond the average position. These locations 
cannot be explained by a comparably distant magnetopause whose position varies 
only as the one-sixth power of the momentum flux but must rather be 
due to an enhanced standoff distance for the bow shock. These distant 
shocks are found to be weaker than average bow shocks on the basis of a 
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decreased jump in the tangential field strength across the shock. An 
unusually low Alfven Mach number is often observed at these times and is 
consistent with an enhanced standoff distance as proposed by Spreiter and 
Jones (1963)* This result agrees with the conclusions of Heppner et al. 
(1967) but the observations extend the position of distant shacks to far 
greater distances. 
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Figure Gaptions 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Position of the bow shock in the solar ecliptic plane 

as determined by measurements on 5 IMP spacecraft. Crosses 

represent the average location on individual passes and the solid 

line hyperbola represents the best fit curve to the points. 

Points have been rotated by 4° to remove the effects of 

aberration due to the earth 1 s motion about the sun. The 

line segments beyond the average shock position represent 

the positions of unusually distant bow shock locations. 

Position of the magnetopause in the solar ecliptic plane 

as determined by measurements on 6 IMP spacecraft. Crosses 

represent the average location on individual passes and the solid 

line ellipse represents the best fit curve to the points. 

. ° 

Points have been rotated by 4 to remove the effects of 
aberration due to the earth* s motion about the sun. A 
solid line has been joined to the ellipse at X = -15 to 
represent the boundary of the tail. 

Relative occurrence frequency of solar wind Alfven Mach 
number M^ and gasdynamic Mach number M as calculated from 
the IMP 4 interplanetary magnetic field and plasma measurements. 
Position of the magnetopause in the noon meridian plane. 

Crosses represent the observed locations within 15° of 
the noon meridian and curves represent the theoretical 
predictions of Olson (1969). 


Figure 5 


Figure 6 


Figure 7 


Position of the magnetopause in the dawn; and dusk meridian 
planes. Crosses represent the observed locations 
obtained within 15° of the dawn and dusk meridian planes 
and curves represent the theoretical predictions of Olson 

(1969) . 

Subsolar magnetopause and bow shock positions during 1967 . 
The histogram trace represents the positions predicted from 
the interplanetary hourly average measurements and the 
crosses and dots represent the observed average positions 
on each orbit. Observed positions have been normalized 
to the subsolar point by taking variations from the average 
curve and plotting them relative to the average subsolar 
points represented by the horizontal lines. 

Distribution of observed and predicted subsolar magnetopause 
positions. Predicted positions have been adjusted so that 
they are centered on the observed average position au 10 . 9 R 
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